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2° : Abstract 

j-H . Laser experiments with optically excited frozen gases entail the excitation of polar- 

C , ization waves. In a continuum approximation the waves are dispersionless, but their 

c/3 I frequency depends on the direction of the propagation vector. An outline is given of 

' the theory of transient phenomena that involve the excitation of these waves by a 

' resonant dipole-dipole transfer process. 

S: 

T3 ■ 1. INTRODUCTION 

a ■ 
o ■ 

O I If the energy associated with the translational temperature of a gas is low compared to the 

relevant excitation energies and to the energy associated with the inverse observation time, 

CN . then the gas can be considered as "frozen" and is analogous to an amorphous solid, or a 

^ . glass. Such a frozen gas can be produced by laser cooling, and the atoms or molecules can be 

i selectively placed in certain excited states. When only a few states are involved, the system 

\ is analogous to a spin glass and can support polarization waves (dipolar Frenkel excitons) 

O ' that are analogous to spin waves. In this paper, we introduce a basic example of such waves 

^ ■ and discuss their excitation. 
G^ ■ 

. In laser experiments, the frozen gas is set up in a non-equilibrium electronic state, and 

^ ! the evolution of this state is measured after some elapsed time. Relaxation processes can 
G I be neglected for the duration of the experiment. Thus, the electronic state evolves towards 
'1^3 ■ dynamical equilibrium through the excitation of polarization waves in the random medium of 
^ ■ the frozen gas. The theory of this process must determine the spectrum of the polarization 

O . waves as well as their excitation probability during the transient evolution to dynamical 
^ ! equilibrium. The two aspects of the problem will often be intermingled. However, we have 

I chosen to work on a system where they can be separated to some extent. 
^ ■ The system under consideration is fully described in an earlier paper by J. S. Frasier, V. 

■ — ' Celli, and T. Blum,0 henceforth to be referred to as FCB. In the present paper, we develop 
an improved theory of the particular case that is discussed in Sec. Ill of FCB, under the 
heading "Sparse ss' — > pp' in a bath of sp — > ps." Simply put, we consider one atom initially 
in state s' that is immersed in a frozen gas of atoms in state s. The s' state is connected to 
a p' state by a dipole matrix element fj,', and an s state is connected to a p state by a dipole 
matrix element /i. The p and p' states are singled out from all possible states, because the 
ss' — >• pp' transition is resonant, i.e., e'^ + es = e'p + e^, and it is known that the transition 
rate peaks sharply at resonance.lii Deviations from resonance are not considered in this 
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paper. The resonance condition is achieved by the apphcation of a static electric field. The 
ss' pp' process excites polarization waves described by sp — ps transitions in the gas 
of s atoms. Typically, the number of s atoms in the volume Q is such that N/Q ~ 10^ 
atoms/cm^. The s' atoms are "sparse" if their number A^' is such that N'/N <^ 1, and in 
that case the result for A^' atoms can be obtained by computing the result for one s' atom, 
averaging over distributions of the s atoms, and then multiplying this result by N'. 



2. THE EQUATIONS OF TIME EVOLUTION 

We begin with the set of coupled differential equations in Eq. (24) of FCB, Fourier trans- 
formed in time: 

^ao = "^^€1 + i, (la) 

I 

UCl = Viao + UlmCm- (lb) 

m 

As in FCB, {t) is the amplitude of the state in which the atom at the origin is in state 
s' and all other atoms are in state s, while q {t) is the amplitude of the state in which the 
atom at the origin is in state p', the atom at position is in state p, and all other atoms 
are in state s. The quantity Vi is the interaction potential 



r (r) = ^ (l - 3 cos^ 9r) = -2^P2 (cos 9^) , (2) 
evaluated at r = r/, and 9r is the polar angle. Similarly, Uim is 

2 2 

f/ (r) = ^ (l - 3 cos^ 9r) = -2^P2 (cos 9^) , (3) 



evaluated at r = r; — rm- The angular factor present in V and U has been chosen for 
mathematical convenience: the actual angular factor is more complicated, as discussed in 
Sec. In FCB, all the results depend only on an angular average of V'^ and, as is shown 
in the Appendix of that paper, one can work with the effective angle-averaged potential 
V (r) = {A/Sy/Sj /i/i'/r^. In the treatment presented here, however, the angular dependence 
plays a more important role. 

The potentials of Eqs. (^ and (^ have particularly simple Fourier transforms. Defining 

V{k) = f (fr V (r) e-'^ \ (4) 



we find, using the standard expansion of e in spherical Bessel functions and spherical 
harmonics,EI 

°° dr 



V (k) = -87r/i/i'P2 (cos 9k) / —32 (kr 

Jo r 



~Hn'P2{cos9k) . (5) 
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Similarly, 



f/(k) = -^f,^P, (cose,). (6) 

With these expressions we can obtain explicit solutions for the polarization waves in the 
continuum approximation, and we can then use these approximate solutions to describe the 
dynamics of a frozen gas state. 



3. THE CONTINUUM MODEL 

We change the sums in Eqs. (1) to integrations over space, so that the evolution equations 
become 



iV f 

^^0= -Q J ^ ^ + ^' ( 

r 

ooc (r) = V{r)ao + -^J d\' U {v - v') c (r') , (7b) 
where c (r) is now an averaged quantity. Fourier transforming, we obtain simply 
t 

VI J (27r) 



^«o = -y 77-3^(k)c(k)+z (8a) 



a;c(k) =\/(k)ao + -f/(k)c(k). (8b) 

The last equation shows that the frequencies of the polarization waves are 

uj (k) = (k) , (9) 

independent of the magnitude of k, but dependent on the direction, according to Eq. (^. 
This peculiar dispersion relation, which is analogous to well-known results in spin wave 
theory, is a consequence of the long range nature of the dipole-dipole interaction. We expect 
it to be useful for values of k up to the order of the inverse of the interatomic spacing.! The 
solution of Eqs. (8) gives 

"°^"^ = a;-M(a;) ' ^'^^ 
where 

, N f (fk Vikf ,,,, 
Miuj) = — o ^4—. 11) 

Carrying out the angular integration, we find 



M(.;) = -A(^f 
6n 



''^ dkP 







m(z), (12) 
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where, leaving for Sec. || a discussion of the effective cutoff kc 



and z = uj/uq with 

no = y/.^^. (14) 
A plot of Re m and Im m on the real axis is shown in Fig. |l[ For large z, 

As expected, m (z) is an analytic function in the upper z-plane, with a cut extending 
from z = — 1/2 to 2 = 1. On the real axis outside the cut m{z) is real, and the spectral 
density on the cut is given by 

vr z"^ 

\m[m{z)] = -y= (16) 

V6 Jz + i 



4. CONTINUUM APPROXIMATION FOR THE RANDOM MEDIUM 



To complete the evaluation of the self-energy M (u) of Eq. (pT^), it is necessary to introduce 
a cutoff in the k integral. We can do this in a way that makes the result applicable, at least 
approximately, to a random medium. 

We know from Sec. II of FCB that the solution of Eqs. (1) for f/ = is of the form of 
Eq. (|10D, with 



Miuj] 



u=o 



(17) 



while in the continuum limit, using Eq. 
N f d^k Vi}if 



u=o 



ijj 



32 iV 
45 



dk ¥ 



(/i/i') 



Therefore, we see that if we make the substitution 
§/>^^^-fE^i^.(cos^.)f, 



(18) 



(19) 



then in the limit of U going to zero our continuum result will agree with the result of FCB. 
Applying this substitution, we find that Eq. (|T2p becomes 
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5m (z 



, (20) 

We can now average over atom positions by the technique described in FCB, Eq. (9), 
with the resuh that, for Imz > 0, 



{(^o)av = L d(3 exp ( i(3uj - (^) 



dp 

'0 Uo V Mo 

where, as in FCB Eq. (A2) 
167r3/2 



exp (iPz - —^-5ipm{z)] , (21) 



/i/x', (22) 



so that 

^ = (23) 

Mo 9 /i 

The integral in Eq. ( pi]) can be expressed in terms of the complementary error function.0 

To develop approximations that automatically include a momentum cutoff we introduce 
the Green's function of Eq. (lb), which satisfies the equation 



Gin — l^in + X! ^ImGmn- (24) 



Then the solution of Eqs. (1) is of the form of Eq. (0), with the exact self-energy given by 
M{u) = Y^yiGimV^- (25) 



Irri 

Similarly, the on-site Green's function, which clearly cannot be treated in the continuum 
approximation, can be written as 

Gnn = fj 7^ (26) 

^ 2—imp^nm^mp,{n]'~'pn 

while, for I ^ n, 

Gin = Gii^[n]UlnGnn + ^ Gim,[n]UmnGnn, (27) 

where G/m,!™] is the Green's function for a medium where the n atom is absent. These 
equations are exact, but we now replace all occurrences of Gnn on the right hand side by 
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(G), where (G) denotes the average value of G„„, which does not depend on n. Averaging 
and Fourier transforming the last equation and setting 

/d h 
— -3G(k)exp[zk-(r,-rJ] (28) 

we obtain 

G(k) - ^^)'^(»^) (29) 
l-{N/n){G)U{k) ^^^^ 

= ^ ( ^ - (G)) 

N \{G)-^ - {N/n)uik) ^ 7 ■ 

Similarly, writing the approximate version of Eq. ( PBj ) in the form 
1 



d(3 exp -i^Yl UL (G) ) (30) 



X exp I iP UnmGmp,[n\Upn | ; 

m^p 



and expanding the second exponential, we obtain the self-consistent equation 
(G) = - r d(3 exp (i(3uo + {-i(3 {G)f'^ u) (31) 
X (^l-2./3/^f/e,,(k)^G(k) + ...), 

where f/g// (k) is the Fourier transform of 

f/e//(r) = t/(r) exp i-if} (G) U\r)) . (32) 



This shows how a cutoff is automatically introduced into the theory. With these values of 
(G) and (Gim) we have approximately 

Miu) = J2 y' (G) + E ^/ (Gir,) K. (33) 

I ly^n 

The self-consistency equation ( ^2|) has been discussed extensively in the approximation 
where the contribution of G (k) is neglected (i.e., where contents of the last set of parentheses 
in Eq. (|32[) are replaced by unity) .0 It is believed that this approximation is adequate when 
the polarization waves are localized in the random medium. 

If (G) is known as a function of u, one can obtain (ao) directly. Comparing Eq. (pHf ) 
with 

we see that for an infinite medium 
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(35) 



It is also possible, in principle, to compute the average of |ao(t)| for the frozen gas. 
The result for ao (uj) with the M (u) of Eq. ( pOD can be inverse Fourier transformed and 
squared, and then the result must be averaged over all r. We have already discussed the 
cut of m (z) in Sec. 0, and one can see from Eq. (p!OD and the plots of Fig. |I] that gq (t) 
will have a pole at some z < —1/2, a pole at z = 0, and a pole at some z > 1. Therefore 
the inverse Fourier transform can be computed by summing the residues of these poles and 
integrating around the cut, a process which is amenable to numerical techniques. However, 
the result is unreliable because the pole at z = is not really present in a random medium. 
(In the continuum approximation, waves of zero frequency travel in the direction where 
P2{cos6) = and are not excited by the coupling V, which vanishes at this angle.) A 
similar calculation using the result of Eq. (|33| ) is very cumbersome. In practice, we have 
found it preferable to solve the equations (1) numerically for a given configuration of a few 
hundred atoms, and then average over an ensemble of configurations. These results will be 
reported in a separate publication. 



5. DISCUSSION AND CONCLUSIONS 

In the course of our theoretical investigation of the transient behavior of an optically excited 
gas, we have come across a particular class of polarization waves. They are Frenkel dipolar 
excitons, with the difference that the degeneracies of the s and p states have been lifted by the 
application of a static magnetic field. First, we note that the spin-orbit coupling separates 
the pi/2 states from the P3/2 states. Then, considering the pi/2 states, we see that the 
interaction between the |s, ]) and pi/2, mj = states is indeed proportional to P2 (cosO), 
where 6 is the angle between the interatomic vector and the applied magnetic field. In the 
actual experiments,!'! there is a spatially varying magnetic field in the magneto-optical trap, 
and the full analysis is much more complicated. 

In the simplest case these waves are dispersionless and their frequency, given in Eqs. (^ 
and (^, depends on the angle 9k between the propagation vector and the field's direction 
according to 1 — 3cos^6'fc = — 2P2 (cos6'fc). These waves have not (to our knowledge) been 
previously recognized, perhaps because their occurrence is confined to systems that can be 
described by the highly simplified Eq. (la), which additionally is treated in the continuum 
approximation. 

We remark that the continuum approximation becomes exact in the long wavelength 
limit and is expected to have a wide range of validity because the dipole-dipole interaction 
is long-ranged. The polarization waves also exist on a regular lattice, with a dispersion 
relation that can be calculated.i'0 In a frozen gas the excitation spectrum will spread out, 
but localization effects are probably small, again because of the long-range nature of the 
dipole-dipole interaction. 

Although in our model system almost all atoms are in an excited s state, they could just 
as well be in the ground state. What is important to the model is that this state couples 
effectively to a single p state. This is a good approximation for the Rydberg atoms of the 
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system we are modeling,Bcl but it is not a good approximation in general. Also, the size of 
the dipole matrix element /i is particularly large for Rydberg atoms. 

Finally, we mention that the dipole-dipole interaction does not include retardation 
effects. However, these are expected to influence the results only by inducing some dispersion 
at very long wavelengths. 

In this paper we consider a particular mechanism of excitation of the polarization waves, 
namely, we allow the resonant transfer ss' — *■ pp' to create the initial p state that then 
interacts with the s states.! However, other mechanisms of creation of the p state would 
work just as well as long as a single p state is created. We have also assumed that if only 
a few p atoms are created, then each can be treated independently in accordance with the 
sparse hmit. 

We outline (but do not carry out in detail in this paper) the calculation of the transient 
behavior involving the excitation of the newly-found polarization waves in a frozen gas. 
While this allows us to establish contact with existing or currently planned experiments, 
one should perhaps think instead of more direct ways of observing the polarization waves. 
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Fig. 1. The real part (solid line) and imaginary part (dashed line) of the function m{z) for z 
just above the real axis. 
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